Electron-phonon interaction in the normal and superconducting states of MgB2 
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For the 40K-superconductor MgB2 we have calculated the electronic and phononic structures 
and the electron-phonon interaction throughout the Brillouin zone ab initio. In contrast to the 
isoelectronic graphite, MgB2 has holes in the bonding cr-bands, which contribute 42 per cent to 
the density of states: A'' (0) — 0.355 states/(MgB2 -eV-spin). The total interaction strength, A = 
0.87 and Xtr ~ 0.60, is dominated by the coupling of the a-holes to the bond-stretching optical 
phonons with wavenumbers in a narrow range around 590 crn^^. Like the holes, these phonons 
are quasi two-dimensional and have wave- vectors close to FA, where their symmetry is E. The 
TT-electrons contribute merely 0.25 to A and to Xtr- With Eliashberg theory we evaluate the normal- 
state resistivity, the density of states in the superconductor, and the B-isotope effect on Tc and 
Ao, and find excellent agreement with experiments, when available. Tc=40K is reproduced with 
fi* — 0.10 and 2Ao/ ksTc = 3.9. MgB2 thus seems to be an intermediate-coupling e-ph pairing 
s-wave superconductor. 



The recent discoveryS of superconductivity with 
Tc^39 K in the graphite-like compound MgB2 has caup 
hectic activity. Density-functional (LDA) calculations 
show that, in contrast to intercalated graphite (Tc < 5 K) 
and alkali-doped fullerides, A3C60 (Tc < 40K) , in MgB2 
there are holes at the top of the B-B bonding a- 
bands, and that these couple rather strongly to opfti-, 
cal B-B bond-stretching modescl with wavenumbersB'cl 
around 600 cm~^. These are the same type of modes 
as those believed to couple most strongly to the tt- 
electrons in graphite and CgoriH'liere their wavenumbers 
are 2.5 times larger, howeverErLl. Rough estimate^lj of 
the electron-phonon coupling-strength for s-wave pair- 
ing in MgB2p,yield: A ^1. Measurements— pL-±he B- 
isotope effectO on Tc, tapaellingO, transportlljil^, ther- 
modynanain^ropertieaijIiil'Ej, and the phonon density 
of statesE^lHl confirm that MgB2 is most likely an 
electron-phop.aa_mediated &avave superconductor with 
intermediateBollZl or strongtj coupling. 

In order to advance, detailed comparisons between ac- 
curate results of Eliashberg theory and experiments are 
needed. Consider again the example of A3C60, also be- 
lieved to be conventional s-wave superconductorsQ with 
Tc's described by the McMillan expression: 
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The LDA values for A are 0.4-0.6 and the value of the 
Coulomb pseudopotential fi* to be used in (^ is pre- 
sumably considerably— larger than the usual value, 0.1- 
0.2, for sp-materialaia due to the small width (~0.5 
eV) of the tt ii,j-subband compared with the on-ball 
Coulomb repulsionQ. For MgB2, where the 7r-band is 
15 eii, broad, one expects fi* to be 0.1-Q-]2 and the 
LDAo plus generalized-gradient correctionE3 to give A 
with an accuracy better than 0.1. The result of such a A- 



calculation will be presented here and should allow us to 
reach conclusions about the superconductivity in MgB2. 

MgB2 consists of graphite-like B2-layers stacked on- 
top with Mg in-between. The primitive translations 
are: a=(V3/2, 1/2, O) a, b=(0,a,0), c=(0, 0, c) , withO 
a — 3.083A= 5.826 ao and c/a — 1.142. In reciprocal 
space, and in units of l-nja, the primitive translations 
are: A=(2/a/3, 0, 0) , B=(-1/a/3, 1, 0) , C=(0,0, a/c), 
and the points of high symmetry are: r=(0,0, 0), 
A=(0,0,a/2c), M=(l/y3,0,0) , L=(l/y3, 0, a/2c) , 
K=(l/v^, 1/3,0) , H= (1/^3, 1/3, a/2c) . To reach a nu- 
merical accuracy exceeding 0.1 for A requires careful 
sampling throughout the Brillouin-zone for electrons as 
well as for phononSr-jdue to the small size of the cylip-, 
drical cr-hole sheetscl. We therefore used Savrasov'ala 
linear-responce full-potential LMTO density-functional 
method, proven to describe the superconducting and 
transport properties oie.g. Al and Pb with high accuracy. 
The Brillouin-zone integratians were performed with the 
full-cell tetrahedron mcthodE3 with the k-points placed 
on the (A, B,C) /24 sublattice. For the valence bands, 
a triple-kappa spd LMTO basis set was employed and 
the Mg 2p-semicore states were treated as valence states 
in a separate energy window. The charge densities and 
potentials were represented by spherical harmonics with 
I < 8 inside the non-overlapping MT spheres, and by 
plane waves with energies < 201 Ry in the interstitial re- 
gion. 

The resulting electronic structure is uraytically iden- 
tical with that of previous calculationaJu. Near and 
below the Fermi level there are two B pz 7r-bands and 
three quasi-2D B-B bonding cr-bands. The a and tt bands 
do not hybridize when kz=0 and tt/c. The 7r-bands lie 
lower with respect to the cr-bands than in graphite and 
have more fcz-dispersion due to the influence of Mg, the 
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on-top stacking, and the smaller c/a-ratio. This causes 
the presence of ^^(=0.056 light and Po-/i=0.117 heavy 
holes near the doubly-degenerate top along FA of the 
CT-bands. For the density of states at ej? = 0, we find: 
N (0) = N^i iO) + N„h iO) + N^ (0) = 0.048-^0.102-^0.205 
= 0.355 states/(MgB2-eV-spin). 

The a and 7r-bands may be understood and described 
with reasonable accuracy near ep using the orthogo- 
nal tight-binding approximation with respectively the 
B pz orbitals and the B-B two-center bond-orbitals 
formed from the B sp^ hybrids: With two Pz or- 
bitals per cell and hopping between nearest neigh- 
bors only (e^=0.a4eV, i^=0.92eV, and t^=1.60eV), 
the TT-bands ara3: (k) = Sz + 2t^ cos ckz ± 

tz ^Jl -f4cos(afcy/2) [cos(afcy/2) -f cos {akxV3/2)]. The 
bonding cr-band Hamiltonian is: 

Ha (k) = tsp2 — 2t^ cos ckz 

{0 COS7 + rcos (a -I- /3) cos a -I- r cos (/? -I- 7) 1 
c.c. cos/3 -|- r cos (a — 7) > 

c.c. c.c. J 

in the representation of the three bond-orbitals per 
cell, with tspi being the energy of the two-center bond, 
and the integrals for hopping between nearest and 2nd- 
nearest bond orbitals in the same layer being respectively 
ib=5.69eV and t'^ = rif,=0.91eV, and with t^=0.094eV 
being an order of magnitude smaller than t^. Moreover, 
a = ik • a, /? = 5k • b, and 7 = ^k- (b — a) . Along FA, 
a=/3=7=0 so that there is a singly-degenerate band of 
symmetry A with dispersion tgp2 — 2t^ cos ckz — Aitb + 
and a doubly-degenerate band of symmetry E with dis- 
persion tsp2 — 2t^ cos ckz + 2{tb+t'f^). The E-band is 
slightly above the Fermi level and its eigenvectors are 
given in the two inserts at the bottom of Fig. |^. The 
Fermi-surface sheets are warped cylinderso which may be 
described by expanding the two upper bands of H^r (k) 
to lowest order in k"^ + ky = k'^y This yields: 

£o-„ (k) = £0 - coscfe - fcjj/mo-n , (2) 

where £0 = ^sp^ -\- 2 {tb -\- t'f^) = 0.58 eV is the average 
energy along FA and the units of fc^/mo-n and k\\ are re- 
spectively Ry and a^"^ . The light and heavy- hole masses 
are respectively rrici — 4/ {tbO^) = 0.28 and niah — 
4/ (3t|,a^) = 0.59 relatively to that of a free electron. 
For energies so closely below e (F) that Eq. (|) holds, 
the number-of-states function is / fc^^^ (e, kz) dkz/kg^ ~ 
fcf|„ (£,7r/2c) /fcl^ = (£0 - £) m^„/fc|2, and its energy 
derivative is therefore constant: N^-n (s) — "m-an/k^^. 
Here, nkg ^ — {2TTao / a)'^ 2/a/3 is the area, and ksz the av- 
erage radius of the Brillouin zone. Not e that th e cr-sheets 
are quite narrow: kp\\ i (7r/2c) /fcBZ=-\/0. 056/2=0. 17 and 
kpw h (7r/2c) /fcBZ=V0-117/2=0.24. 

The dynamical matrix was calculated for q-points 
on the (A, B,C) /6 sublattice using the B-mass 10.811 
which corresponds to the natural mix of isotopes. The 



phonon dispersions (q) and density of states F {lo) are 
shown in Fig.|l|. The agreement between our F (tbiJ|-and 
those obtained from inelastic neutron scatteringla'tH is 
excellent; our peaks at 260 and 730 cm~^ (32 and 90 meV) 
are seen in the experiments at 32 and 88 meV. For the fre- 
quencies of the optical F-modes we get: 335 cm^^ (Ei«), 
401cm-i (A2„), 585cm-i (Ejg), and 692 cm^^ (Bi^). 
These values are in good agreement with those of previ- 
ous calculations, except for the all-important Y^2g modes 
where the LAPWH and pseudo-potentiaH calculations 
give values which are respectively ^lOOcm^^ smaller 
and larger than ours. These doubly degenerate modes 
are the optical B-B bond-stretching modes (obs). Close 
to FA, they have exactly the same symmetry and simi- 
lar dispersions as the light and heavy cr-holes, although 
with the opposite signs. The E eigenfunctions shown at 
the bottom of Figj^ now refer to displacement patterns, 
e.g., {—1,0,1} has one bond shortened, another bond 
stretched by the same amount, and the third bond un- 
changed. These E displacement patterns will obviously 
modulate the electronic bond energy, tsp2 , in such a way 
that the light holes couple to one, and the heavy holes 
to the other mode, with the same diagonal {gambs) and 
no off-diagonal matrix element. An and PickettJj judged 
that this electron-phonon (e-ph) matrix element will be 
the dominating one. This, we shall confirm, although we 
dispute their estimate of A. 

We now turn to the superconducting properties. Since 
the Fermi surface has two a and two 7r-sheets, one might 
expect anisotropic pairing with d iffer ent gaps on differ- 
ent sheets. The experimental dataOO, however, demon- 
strate that a change of residual resistivity corresponding 
to a change of impurity scattering by two orders of mag- 
nitude hardly affects Tc, although Tc for anisotropic pair- 
ing is very sensitive to impurity scattering. The relative 
linewidth of the mq-phonon due to e-ph coupling isE3: 

^"^ = 27r (5 [£„ (k)] 5 [sn' (k + q)] |g„k,«'k+q,m|^ 

nn'k 

= TTiV (0) uj,n (q) Ki (q) , 

where the factor 2 is from spin degeneracy and 
is the average over the Brillouin zone, so that 
^ (0)=Snk [^n (^)] ■ have (safely) assumed that 
(q) < q • v„ (k) , where v„ (k) = Vk^n (k) is 
the electron velocity. The c-ph matrix element is: 
.9nk,n' k+q,m = ("k \5V\ n'k -|- q) /(5Qmq, where the dis- 
placement in the i-direction of the jth atom is related to 
the phonon eigenvector Cij^mq and displacement JQmq 
by: 5Rij = /^2MjW„q. The Ehashberg 

spectral function is: 

(c) F (c) ^ J2 - c.„, (q)] , 

27rA^ (0) ^ uj,n (q) 

^ ^ mq ^ ' 

and the strength of the e-ph interaction is finally: A = 
2 cj-W {cj) F {cj) dcj = A„, (q) . 
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FIG. 1. Left: Calculated phonon dispersion curves in MgB2. The area of a circle is proportional to the mode-A. The insets at 
the bottom show the two FA E eigenvectors (un-normalized), which apply to the holes at the top of the u-bands (bond-orbital 
coefficients) as well as to the optical bond-stretching phonons (relative change of bond lengths). Right: F [uj) (full curve and 
bottom scale), {uj) F {ui) (broken) and q^^ (oj) F (cj) (dotted). See text. 



The dominance of the a-<j coupling via the optical 
bond-stretching mode is clearly seen in Fig. |l|, where 
the area of a black circle is proportional to Am (q) ■ 
Along FA, except when q • v < oj, only the small k^- 
dispersion (f^) makes \m (q) not diverge so that the nu- 
merical values are inaccurate due to the relative coarse- 
ness our k-mesh. The nearly cylindrical cr-sheets, whose 
diameters are of about the same size as the smallest, 
non-zero q\\ on the affordable (A, B, C) /6-mesh, re- 
quire even more care in the numerical q-integration: In 
case of a single cylindrical sheet with p holes, A((7|[) 
has the well-known Imx (g||,Lj ^ 0)-form: \{q\\^/\ = 

(2ttpx^/1 — x'^) 9 {1 — x) with x=q\\/2kp\\. This func- 
tion vanishes when gn >2fc^|| , has a flat minimum of value 

(irp) ^ near q^\=y/2kp^^, and has integrable divergencies 
at (7||=2fc^|| and 0. The proper average of A is A. This 
means, that A (q) calculated on a coarse mesh scatter vi- 
olently for small |q| , but that weighting with A/A(<7||) 
gives the same, correct result for all these points, pro- 
vided that warping, as well as k, k'-dependence of g and 
Lu, are neglected. In case of two cylindrical sheets, and no 
coupling between them, A„ (g||) /A„ should be weighted 
by m^/ (rnf + m^) . In our numerical evaluatimi of the 
e-ph interaction with the linear-response codell3 we dis- 
carded the values of A„j (q) with q along FA, and 
added those on the (A/12, B/12, C/6)-mesh for which 
V^kpin < < V2kpiih. The result was: A=0.62-h0.25, 
where 0.62 was the contribution from q's so small that 
a-a coupling occurs, and 0.25 was the contribution from 



the remaining part of q-space, which must involve a tt- 
sheet. Had we included the inaccurate Am (q)-values 
along the FA-line, the a-cr result would have been 0.72 
instead of 0.62. The result was finally checked by us- 
ing the approximate A(g||) /A correction for the point 
q=A/12. This yielded 0.58 instead of 0.62. In conclu- 
sion: A = 0.87 ± 0.05 = (0.62 ± 0.05) + 0.25 = A^ -I- A^. 

The Eliashberg function shown on the right-hand side 
of Fig. ^ is dominated by the large a-a peak around 
^ofcs=590 cm^^=73 meV. The facts that the tr-sheets are 
narrow, warped cylinders whose coupling is dominated 
by intra-sheet coupling via the optical bond-stretching 
mode, and that the coupling between a- and 7r-sheets is 
negligible, lead to the following approximation: 

(lo) F (lo) « al {to) F (lo) [N^ (0) /N (0)] -I- 
\g,^obs\' S {cu - LOobs) [Nil (0) + Nl^, (0)] /N (0) , 

where (w) F (w) is the usual expression, but with tt- 
electrons only. In An's and Picket's estimateO, Acr=0.95, 
a factor [N^i (0) -t- TV^^ (0)] /N„ (0) N (0) = 0.24 appears 
to be missing. The rigid-atomic-sphere estimate A=0.7 
by Kortus et al. is closer to our value 0.87. 

Knowing {lo) F (w) and a value of the Coulomb 
pseudopotcntial /i* (wc), we soljie the Eliashberg equa- 
tion on the real frequency axiaHj, and obtain Tc=40K 
if jj* (a;c)=0.14. Taking retardation effects into ac- 
count, we find /i* = p* [ujc) / [1 + /i* {ujc) In (wc/win)] — 
0.10, where a;in=504cni^^=62meV is obtained from: 
0= In (w/win) W^a^ {u) F {u) duj, and the cut-off fre- 
quency is taken as Li;c=10maxLj=8000cni^^. This value 
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of /X* is at, the lower end of what is found for simple 
sp-metalalj. The relation back to a screened Coulomb 
interaction U is: /x* = /x/ [1 + In (tjp/tJin)] , where 
fi—UN{0) and ujp ~7eV is the plasma frequency given 
below. We thus find: fi— 0.19 and U—l.leV, which are 
normal values. Had we used the approximate McMillan 
expression the slightly higher value /i*=0.14 would 
be needed to reproduce the experimental Tc. 
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FIG. 2. Normalized density of states (full) and the neg- 
ative of its energy-derivative (dotted) as obtained from the 
Eliashberg equation with jj,*— 0.10 and T=3K. 

In Fig|2| we show our Eliashberg calculation with 
^*=0.10 of the density of states, Ns{e)/N{0) = 

Re ' ' 



in the superconductor. The BCS 

singularity is at e= A3K (0)=6.8 meV, which is in accord 
with the 4.9-6.9 meV found in tunneling experimentsO. 
This yields: 2Ao/A:bTc=3.9 which is slightly higher than 
the BCS value of 3.52. The distinct feature near 80meV 
corresponds to the peak in a'^{u!)F{u!) at 73meV, shifted 
by the 6.8 meV gap. The latter function is also shown 
in the figure together with the measurable quantity 
~d^I/dV^ ~ -dNs{e)/de. 

We have calculated the change in Tc upon isotope sub- 
stitution of for and get: 5Tc=~1.7K, which 
corresponds to the exponent — (51nXc/(51nAfB=0.46. 
This agrees well with the measuredB value: STc=- 
1 K. For the change of the gap, which may be mea- 
sured in tunnelling and optical experiments, we calcu- 
late: (5Ao=-1.9cm~^, which corresponds to the exponent 
-(51nAo/51nMB=0.38. 

Finally, we have considered transport properties in the 
normal state. Here, solution of the kinetic equation leads 
to the transport e-ph spectral function a^^ ^{uj)F{uj), 
and similarly for y and z. These components are 
given by the previous expressions, but with the ad- 
ditional factor [vl^ (k) - Vnx (k) Vn'x (k + q)] / (v^) in- 



(k)5[(e„(k)]. In 



serted. (v^) = y^j 

Fig. |l| the directional average, aj^{uj)F{uj), is seen 
to have the same shape as a^(w)F(w), except for the 
a-a interaction via the optical bond-stretching modes, 
whose aj^(w)F(w) is smaller, presumably due to the 
near two-dimensionality of the cr-bands. As a re- 



sult. At,- = 0.60. For the plasma frequencies, oj^ ^ = 
47re2A^ (0) (w^) / [abc] , we find: uj p =tOp,y =7. 02 eV 
and ijjp,z=6.68eV. Also the temperature dependence 
of the specific dc-resistivity calculated with the 
standard Bloch-Griineisen expression, Pdcx{T) = 
{tt/loI^T) J^ujsmh-^ {L0/2T) al {u;)F{u)du, is nearly 
isotropic and, as shown in Fig.H, is-, in accord with 
recent measurements on dense wiresEHI over the entire 
temperature range. The crossover from power-law to 
linear temperature dependence is seen_.to occur near 
maxa;/5=160cm~^=230K, as expectedEJ. 
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FIG. 3. Calculated dc-resistivities in different directions 
compared with the experiment in Ref. 

In conclusion, we have presented an accurate ah ini- 
tio calculation of the c-ph interaction in MgB2 and find 
A = 0.87 ± 0.05. Eliashberg theory with /i*=0.10 gives 
good agreement with available experiments and several 
predictions. The unexpected high Tc is due to the large 
A- value caused by the presence of holes in the B-B bond- 
ing cr-band and the relative softness of the optical bond- 
stretching modes. MgB2 thus seems to be a simple and 
clear case of an intermediate-coupling e-ph pairing s-wave 
superconductor. 
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